Magnetic stress as a driving force of structural distortions: the case of CrN 
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We show that the observed transition from rocksalt to orthorhombic P„ma symmetry in CrN can 
be understood in terms of stress anisotropy. Using local spin density functional theory, we find that 
the imbalance between stress stored in spin-paired and spin-unpaired Cr nearest neighbors causes 
the rocksalt structure to be unstable against distortions and justifies the observed antiferromagnetic 
ordering. This stress has a purely magnetic origin, and may be important in any system where the 
couphng between spin ordering and structure is strong. 
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The simultaneous occurence of structural distortion 
and change of magnetic ordering is observed in several 
magnetic compounds. In general, the interplay between 
structural and spin degrees of freedom depends on the de- 
tailed configuration of the electronic and phonon states, 
and a unique picture that elucidates the coupling mech- 
anism is missing. A large class of systems where this 
interplay is strong are materials containing ions with de- 
generate (usually eg) orbitals (Jahn- Teller ions), such as 
Mn'^+, Cr^+ or Cu^"*" The degeneracy causes insta- 

bilities that can be relieved by structural distortions or 
orbital (and eventually charge) ordering. In this paper 
we focus on CrN, one of the most widely employed mate- 
rials in coating technology, due to the extreme hardness 
of CrN films. It is also the prototype material of strong 
magneto-structural interactions, since here, more than in 
any other case, the coupling is clearly revealed: at room 
temperature it is paramagnetic (PM) in the RS struc- 
ture, like most of the transition metal nitrides, but at 
Neel temperature rAr=273-286 K, a both structural and 
magnetic transition to an antiferromagnetic (AFM) or- 
thorhombic P„ma phase is observed ||-||]. 

But other compounds presents indications of strong 
spin and structural coupling. To name a few of them, 
in the perovskite LaMnOa the distortions are necessary 
to stabilize the observed A-type antiferromagnetic order- 
ing. The doped cobalt perovskite Lai-ajSr^jCoOa exhibits 
unusual properties as a function of doping and temper- 
ature 01, and in the diluted magnetic semiconductors 
Zni_a;Crj,Se, the magnetic Cr^'^ impurities are strongly 
affected by Jahn- Teller distortions Also large is the 
class of materials manifesting magnetoelastic effects, such 
as changes of magnetic anisotropy coupled to distortions 
and strains. For example, the martensitic phase transi- 
tion of Ni2MnGa from cubic Heusler to tetragonal struc- 
ture is accompanied by a change in magnetic anisotropy 
and driven by magnetoelastic coupling. |^ . In the layered 
Cu/Co/Ni/Cu/Si(001) the magnetic anisotropy can be 
induced by strain and controlled through the film thick- 
ness [0, and in CeNiSn the magnetic ordering can be 
changed by application of uniaxial strain JTl]] . 

For all of these materials, the stress is an important 



quantity, being naturally able to bridge electronic and 
structural properties. The stress relief argument is very 
popular in surface physics to explain processes like sur- 
face reconstructions, adsorbtions and growth: after a 
structural perturbation, a stress excess is produced in the 
system (e.g. the surface stress), that is relieved by further 
structural rearrangements involving symmetry breaking 
and eventually change of atomic density at surface. In 
this paper we show that is possible to apply similar ar- 
guments to magnetic systems. Specifically, we use the 
concept of stress relief to explain the stability of one mag- 
netic ordering with respect to another, and to show how 
the change of ordering is related to a structural distor- 
tion. 

In analogy to the surface stress, we define magnetic 
stress to be the stress associated with a change of spin 
ordering. If and represent two magnetic phases of 
a system with the same structure, the magnetic stress T 
associated with the transition is: 
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where e is the strain tensor. T has a purely magnetic 
origin, since it originates from changes in spin orienta- 
tions, not structural transitions, but may be relieved by 
structural transitions. We also define the magnetic stress 
per bond by mapping T onto the Heisenberg model. 
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where tij is the strain derivative of the usual energy 
exchange interaction parameter Jy , and represents the 
change of stress caused by turning the spin orientation 
on atoms i and j from parallel to antiparallel. A posi- 
tive tij means that the single spin-flip produces a tensile 
stress that can be relieved by shortening the i-j distance. 
A negative tij, instead, represents a compressive stress 
that is relieved by stretching the i-j distance. The mean- 
ing of the definitions of Eqs.|l| and |^ will be transparent 
when applied to the practical case of CrN. 

Below Tn, CrN is orthorhombic, with an unusual mag- 
netic ordering |] (see Figj^, left) built up by (iTO) FM 
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layers whose spin direction alternates up and down after 
each 2 layers, so that the system is antiferromagnetic 
(AFM) overall (in the following this phase is labeled 
AFM^^-i^j). Along with the AFM^^-^j, we consider two 

other AFM phases as possible competitors: the AFM^I^^^-^j , 
similar to the former but with a one-by-one compensation 
of FM layers along [iTO] (Fig.|l|, right), and the AFM^^^^j, 
made by (111) FM layers alternated along [111], com- 
mon for some strongly ionic transition-metal oxides (e.g. 
MnO). 

Our calculations are based on local spin density func- 
tional theory (LSD) ||l^, and employ plane- wave ba- 
sis and TrouUier-Martins pseudopotentials |l^ (this 
methodology already described successfully some prop- 
erties of CrN ]l^-p7[). The microscopic stress tensor 
is calculated in LSD ||l^ (within the framework of the 
ABINIT [|l9| project) and problems of finite-size basis 
have been efficiently circumvented pOt . 

To rationalize the RS-to-orthorhombic transition, we 
can picture it as a three step process: first, an hydro- 
static expansion of the volume accompanies the non spin 
polarized-to-spin polarized transition. This is the dom- 
inant step in terms of energy stabilization. Second, a 
further stabilization is obtained from the competition of 
different magnetic orderings with the same (cubic) struc- 
ture. Third, a shear strain is appHed onto the (001) plane, 
in accordance with the experimental finding |^ . The ob- 
served distortion consists of a reduction of the lattice 
parameter a and an increase of b (of ~ 2%) with respect 
to their ideal values (a = ao-\/2, b = bo/^/2, and c = ao, 
where ag is the cubic lattice constant), in a way that the 
Cr-N distance is unchanged. Although a second-order 
effect with respect to the hydrostatic expansion, we will 
show that the distortion is fundamental in stabilizing the 
observed AFM^^-^^j phase. 

The effects related to the first step are described in 
Fig.^, where the energies and stresses of PM and FM or- 
derings in the RS structure are shown. In the PM phase, 
the high density of states (DOS) at the Fermi energy 
(Ef) strongly drives CrN towards a Stoner instability. 
Although the high DOS peak at Ef disappears in the 
FM phase, the spin splitting docs not manage to open 
a gap, thus, FM CrN is still a (poor) metal ||l^. It is 
also a robust magnetic compound, with a magnetic mo- 
ment (M=2.09 hb per Cr atom) almost entirely due to 
the Cr t2g splitting, and a chemical configuration close 
to Cr2+N2-. 

We find a difference of 0.182 eV/formula unit between 
PM and FM energies, calculated at the equilibrium aq of 
the respective phases, that are 7.72 a.u. and 7.84 a.u. for 
PM and FM, respectively. Thus, after polarization, the 
lattice constant increases ~ 1.5%. This strong magne- 
tovolume effect is related to a large magnetic stress. At 
the PM lattice constant (where T^*^ is zero by defini- 
tion), the spin-polarization produces a compressive stress 



rpFM_-[^ gg eV/formula unit, whose relief corresponds to 
an energy gain A£'^*^=T^^^ A a/a~ 0.03 eV/formula 
unit. 

The increase in compressive stress is a consequence of 
some charge localization (mostly Cr d t2g charge) around 
the Cr nucleus, due to the spin polarization. The charge 
localization increases the kinetic stress, and, in turn, de- 
creases the tensile electron-ion interaction. 

Within this expanded volume, but still retaining the 
RS symmetry, we consider the competition among differ- 
ent magnetic orderings (the calculated lattice constant is 
quite similar for FM and AFM, i.e. a change of magnetic 
ordering does not produce that significant hydrostatic 
compression) . We find (Table ^ all the AFM phases are 
favored over the FM one. Moreover, the AFM|!^j^-^j, and 

not the observed AFM^^j^^j , is the stable phase in the RS 
structure. 

The density of states (DOS) of a single Cr ion in the 
AFM^^-^j phase is shown in Fig. ^. The Cr configuration 

is nearly Cr^+, with magnetization M=2.12 /i^. Notice 
that in the (001) plane of the orthorhombic cell, the t2g 
orbitals (whose DOS is drawn by the shaded area in the 
Figure), not the e^s, are directed along x and y. The Cg 
orbitals are non degenerate, due to the AFM symmetry, 
and each of them is occupied by nearly half an electron. 
Their total contribution to the magnetization is only ^ 
0.2 fjLB- The exchange splitting for t2g orbitals is ^ 2 
eV, with the t\g almost completely filled. The global t2g 
occupation is 2.76, and the contribution to the magnetic 
moment is 1.9 ^lb- 

We can understand the competition among spin order- 
ings by a two-parameter Heisenberg model, taking into 
account the nearest neighbor Cr-Cr (Ji), and the next 
nearest neighbor Cr-N-Cr interactions (J2). The calcu- 
lated parameters are Ji=-9.5 meV, and J2=4 meV. For 
the Cr-N-Cr interaction there is competition between su- 
perexchange (AFM) and double-exchange (FM) mecha- 
nisms involving Cr Cg and N p orbitals {dpa hybridiza- 
tion). The resulting coupling is slightly FM (i.e. J2 is 
positive). The direct Cr-Cr coupling, instead, is medi- 
ated by t2g orbital interactions {dda), and is AFM (i.e. 
Ji is negative), as we may argue assuming, as usual, 
the Hund's rule. The formal equivalence of Cr^+ and 
Mn'^"'" ions may suggest some similarities between CrN 
and LaMnOa. Indeed, in both of them the magnetic ion 
is surrounded by an octhaedral cage of cations. However, 
in LaMnOs the Mn coupling via oxygens is the only rele- 
vant interaction, whereas, for CrN the Cr-N-Cr coupling 
is smaller than the direct Cr-Cr AFM coupling. As a con- 
sequence, the stability increases with the number of AFM 
nearest neighbors (8 for the AFM^^-^j , 6 for AFM^^-^j and 

AFM[iii]). Also, due to the symmetry of the Cr fee sub- 
lattice, no other arrangement with more than 8 AFM 
nearest neighbors is possible. In other words, the signs 
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of Ji and J2 rule out the possibility that other orderings, 
not taken into account here, might be more stable. 

On top of the RS structure, we finally apply a shear 
strain e onto the (001) plane of the competing phases 
AFM^^-gj and AFM^^-^^^ (Table ||). First consider the 

stress values for the undistorted (i.e. e = 0) structure: 
the AFM^^-^j phase is under a condition of appreciable 

stress with respect to the AFM^^-^j. This stress orig- 
inates from the spin ordering asymmetry on the (001) 
plane: in the AFM^^-^j phase, each Cr has two spin- 
antiparallel nearest neighbors along x, and two spin- 
parallel along ij. The bonds between AFM neighbors are 
under condition of tensile stress (i.e. T^^^ > 0), whereas 
the bonds between FM neighbors store compressive stress 
(i.e. Tyy < 0). Clearly, a planar distortion will be able to 

relieve the stress and lower the energy of the AFM? - , . 

°-' [110] 

In contrast, the AFM^^-^^j, which is isotropic on (001), 
cannot undergo the same stabilization: each Cr has two 
spin-parallel nearest neighbors in both x and y directions, 
and the eventual energy gain due to the stretch along y 
would be lost due to the compression along x. Also, T^z 
is different for the two phases: in the AFM^^-^^j , the larger 
number of FM Cr-Cr bonds having an orthogonal com- 
ponent causes an excess of compressive stress along z. 

Applying Eq.^ we can express the stress (at e^O) in 
terms of exchange-interaction parameters. For a planar 
strain on (001), we have three independent parameters 
corresponding to the in-plane Cr-Cr bonds (ti), and the 
planar and orthogonal components of the out-of-plane 
Cr-Cr bonds {t2x and t2z)- Our results gives ti=1.04 eV, 
t2x= 0.24 eV, and ^22=0. 20 eV. Thus, the parallel-to- 
antiparallel spin flip causes a significantly large tensile 
stress for the Cr-Cr bonds parallel to (001). 

On the basis of the magnetic stress results at e=0 we 
are able not only to predict the occurence of the dis- 
tortion for AFM^^-^j, but also to estimate the energy 

gain due to the stress relief. Indeed, a 2% distortion 
gives an energy gain AE = {Txx - Tyy) ^ ~ 0.04 eV 
which is enough for the AFM^^-^j to be stabilized over 

the AFMj!^j^-^j . This estimation can be compared with the 

direct calculation of energies and stresses at e 7^ 0. As 
expected, the distortion progressively relieves the stress 
and stabilizes the observed phase. Our results gives a 
theoretical distortion slightly larger than the experimen- 
tal value: at e = 4% most of the planar stress is relieved. 
In contrast, the distortion on the AFM^^^^-^^ increases the 

stress and is energetically unfavorable. 

The stabilization of AFM^^-^j ordering is ultimately 
due to different stress stored in the bonds between 
spin-antiparallel and spin-parallel Cr nearest neighbors. 
Stress asymmetry and ordering of t2g orbitals are inti- 
mately related: since the Cr-Cr interaction occurs via 
direct t2g coupling and is AFM, Ji should vanish in the 



limit of large nearest neighbor distance, whereas a short- 
ening of the distance should increase the overlap between 
t2g orbitals, making Ji even more negative. 

This explains why a bond contraction stabilizes the 
spin-antiparallel interactions, while a bond stretching 
stabilizes the spin-parallel interactions. In Fig.^ we show 
the DOS of ^2-2 and dy-2 orbitals (that contribute most to 
the Cr-Cr planar coupling) for the undistorted AFM^^-^j 

ordering (obviously they are degenerate for AFM^^-^^). 
The tendency of these orbitals to increase the overlap 
(and, thus, the hybridization) along i, and to stay more 
localized along y is clearly visible. 

In conclusion, in this paper we have given a defini- 
tion of magnetic stress, and proposed it as a tool for un- 
derstanding and predicting the subtle relations between 
magnetic ordering and structural properties. For CrN, 
we show that the coupling between the AFM^^-^j order- 
ing and the orthorhombic Pnma structure is driven by 
the relief of tensile stress stored in the AFM Cr nearest 
neighbors, and associated with the direct interaction of 
t2g orbitals. Our results for CrN may stimulate the appli- 
cation of arguments based on stress relief to many other 
magnetic compounds. Indeed, the predictive power of 
the stress should be important in every case where mag- 
netism and structural properties are strongly related. 
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TABLE I. Energies E (in eV/formula unit) of RS CrN for 
different magnetic orderings, calculated at the respective the- 
oretical lattice constants. All the values are relative to that 
of the AFM^^-pj , that is the most stable phase overall, as long 
as the distortion from cubic symmetry is not considered. 
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TABLE IL Energies E (in eV/formula unit) and stresses T 
(eV/formula unit) for AFM^^^-^jj and AFM^^^-^^ CrN as a func- 
tion of the applied (001) planar strain e. The corresponding 
changes of lattice parameters are 5a = -ea, and 5b = eb, i.e. 
the X axis is compressed, and the y is stretched. All values 
are relative to that of the undistorted AFMj^y^j phase. 
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FIG. 2. Energies and stresses (per formula unit) of PM and 
FM phases of RS CrN, as a function of the lattice constant 
ao. 
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FIG. 3. Orbital-resolved DOS of one Cr ion in the AFM^^-^^j 
CrN (undistorted). Due to the 45° rotation of the (001) plane. 
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FIG. 1. Orthorhombic cell of CrN (4 formula unit) in the 
observed AFM^^^j^j spin ordering (left) and in the AFM^^^-^j 
ordering (right). The cartesian axes of the cubic cell (x' and 
y') are rotated 45° with respect to x and y. 
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FIG. 4. 

RS phase. The vertical solid line corresponds to the Fermi 
energy. 
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